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ABSTRACT 


Angle  and  velocity  resolved  inelastic  He  scattering  has  been  used  to 
study  how  the  surface  dynamics  of  thin  rare  gas  films  evolve  on  a  layer-by- 
layer  basis.  Surface  phonon  dispersion  relations  for  ordered  1,  2,  3  and 
25  layer  Kr  films  physisorbed  on  Ag(lll)  are  presented  along  ?M  across  the 
entire  Brillouin  zone.  The  monolayer  data  are  dispersionless,  indicative 
of  an  Einstein  oscillator  mode.  In  comparison,  the  25  layer  film  has  a  well 
developed  Rayleigh  wave,  typical  of  a  thick  crystal  surface.  Excitation 
linewidths  for  monolayer  Kr,  which  vary  across  the  zone,  are  also  briefly 
discussed. 


i 
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The  physical  and  chemical  properties  of  atomically  thin  adsorbed  films 
are  currently  being  examined  in  great  detail  by  a  variety  of  experimental 
and  theoretical  techniques.  In  addition  to  their  technological  importance, 
studies  on  these  systems  are  motivated  by  a  desire  to  understand  the  many 
factors,  including  the  essential  role  played  by  the  underlying  substrate, 
which  contribute  to  the  fascinating  phenomena  associated  with  them,  e.g.  a 
wide  diversity  of  ordered  2-dimensional  structures,  phase  transitions,  and 
novel  catalytic  activity.  To  date,  much  experimental  effort  has  been  directed 
towards  elucidating  the  geometric  and  electronic  properties,  i.e.,  the  static 
properties,  of  these  systems.  Only  recently  have  vibrational  probes,  such 
as  low  energy  electron  energy  loss  spectroscopy  [1],  inelastic  neutron  scat¬ 
tering  [2],  and  inelastic  He  scattering  [3]  allowed  researchers  to  directly 
examine  the  dynamical  properties  of  atomically  thin  films.  Of  these  techniques, 
inelastic  He  scattering  is  the  best  suited  for  examining  low  energy  (<  5  meV) 
collective  excitations  since  it  has  excellent  energy  resolution  (<  250yeV) , 
high  surface  specificity,  and  momentum  sufficient  for  probing  excitations 
across  the  entire  surface  Brillouin  zone  (.SBZ) . 

Tn  this  Letter  we  present  inelastic  He  scattering  data  which  demonstrate 
how  the  surface  phonon  dispersion  relations  for  Kr(lll)  films  physisorbed  on 
Ag(.lll)  evolve  on  a  layer-by-layer  basis,  from  monolayer  through  thick  film, 
across  the  entire  S8Z.  These  model  systems  lend  themselves  to  accurate  theoretical 
analysis  as  a  great  deal  of  information  exists  pertaining  to  rare  gas  inter¬ 
action  potentials  and  rare  gas-metal  physisorption  potentials,  including  two- 
hody  and  higher  order  terms. 

The  dispersion  curves  were  mapped  out  using  angle  resolved  He  time-of- 
flight  (JOF) .  The  measurements  were  performed  in  a  UHV  scattering  apparatus 
(base  pressure  ^2  x  10~^  torr)  equipped  with  a  rotatable  and  differentially 
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pumped  quadrupole  mass  spectrometer  (crystal-to-ionizer  distance  14.45  cm, 
angular  resolution  0.67°).  An  important  feature  of  the  apparatus  is  that 
the  incident  and  final  polar  scattering  angles  are  independently  adjustable. 

Data  can  therefore  be  collected  across  the  entire  SBZ  using  a  fixed  incident 
angle  for  a  given  surface;  in  large  part  eliminating  the  necessity  of  correct¬ 
ing  the  inelastic  transition  probabilities  for  intensity  variations  that  occur 
due  to  selective  adsorption  resonances  in  the  incident  channel.  The  mechanic¬ 
ally  modulated  18  meV  He  beam  (av/v  ^  1%.  shutter  function  %3usec  FWHM)  was 
produced  in  a  liquid  nitrogen  cooled  nozzle  beam  source.  The  Ag  single  crystal 
was  cut  and  polished  to  within  0.5°  of  the  (.111)  plane,  as  confirmed  by  Laue 
X-ray  back-reflection.  Surface  cleanliness  was  checked  using  Auger  spectroscopy. 
When  needed,  the  crystal  was  cleaned  using  Ar+  bombardment,  and  subsequently 
annealed  at  750  K.  Surface  coherence  was  checked  by  monitoring  the  specular 
scattering  lineshape  of  a  64  meV  He  beam,  and  analyzed  using  simple  pertur¬ 
bation  theory  £4],  indicating  a  coherence  length  of  ca.  100  A.  The  azimuthal 
orientation  was  determined  using  H 2  diffraction. 

Our  thin  film  growth  procedures  were  guided  by  extensive  information 
that  was  already  available  about  the  structure  and  thermodynamics  of  rare 
gas  overlayers  on  Ag(lll)  [5].  Ar,  Kr,  and  Xe  form  azimuthally  aligned  but 
transl ational ly  incommensurate  lattices  on  this  surface.  Monolayers  were 
grown  by  controlling  the  pressure  of  a  differentially  pumped  dosing  beam 
and  the  surface  temperature.  Overlayers  were  grown  by  subsequent  computer 
controlled  deposition  on  well  annealed  monolayers.  The  quality  and  struc¬ 
ture  of  the  Kr  surfaces  were  checked  using  He  diffraction,  as  shown  in  Figure 
1  for  four  different  overlayers.  Specular  peak  widths  and  intensities  in¬ 
dicate  that  the  coherence  lengths  of  these  surfaces  are  only  slightly  less 
than  the  Ag  substrate,  i.e,  we  are  not  growing  roughened  films.  The  very 
low  diffuse  scattering  seen  in  Figure  1  also  confirms  the  high  quality  of  the 
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Kr  overlayers.  The  lattice  constant  for  Kr  monolayer,  4.02  A,  is  just 
slightly  larger  than  that  for  the  2,  3,  and  25  layer  films,  4.00  A.  When 
the  experiments  with  an  overlayer  were  completed,  the  coverage  was  confirmed 
using  temperature  programmed  thermal  desorption  (TPD):  a  Kr  tri layer  TPO 
spectrum  is  shown  in  the  inset  of  Figure  3.  Since  the  monolayer  desorption 
peak  occurs  at  a  higher  temperature  than  the  multilayer  peaks,  the  ratio  of 
the  total  integrated  area  to  that  of  the  monolayer  gives  a  Quantitative  mea¬ 
sure  of  the  coverage.  Only  data  collected  on  films  typically  falling  within 
5*  of  a  particular  overlayer  thickness  were  used  in  our  analysis  for  1,  2, 
and  3  layer  structures. 

TOF  spectra  were  taken  using  a  multichannel  scaler  (lusec/chan  dwell 
time),  and  required  approximately  1  hour  of  signal  averaging.  The  first  step 
in  the  analysis  was  to  fit  the  TOF  data  using  a  non-linear  least  squares 
routine.  Flight  times,  referenced  to  the  elastic  specular  TOF,  were  then  used 
to  determine  the  phonon  energies.  The  parallel  momentum  exchange  of  a  given 
peak  could  then  be  calculated  using  the  kinematic  constraints  imposed  by  the 
incident  and  final  scattering  angles,  the  phonon  energy,  and  the  conservation 
of  crystal  momentum.  A  more  detailed  description  of  the  experimental  apparatus 
and  procedures  will  be  given  elsewhere  [6]. 

Figure  2  shows  typical  TOF  spectra  for  the  four  different  Kr  overlayers. 
The  well  defined  peaks  indicate  that  single  phonon  exchanges  are  the  dominant 
inelastic  scattering  channel.  Each  TOF  spectrum  also  includes  an  elastic 
feature,  which  we  attribute  to  incoherent  scattering  from  crystal  defects. 

The  computer  simulations  in  Fig.  2  are  forward  convolutions  which  have  been 
scaled  to  show  excitation  lineshapes,  but  not  probabilities,  and  which  i)  use 
as  input  the  experimental  dispersion  curves  ii)  include  the  instrument  trans¬ 
fer  function  and  iii)  assume  single  phonon  exchanges  having  del ta- function 


energy  widths.  When  the  phonon  excitation  lineshapes  are  allowed  to  have  finite 
energy  widths,  the  simulations  in  most  cases  fit  the  experimental  spectra,  yield¬ 
ing  information  on  phonon  relaxation  [6].  Most  of  the  spectra  were  intentionally 
taken  at  the  same  incident  angle  (9^  *  45°).  An  observation  made  on  Xe  monolayers 
indicates  why  a  constant  incident  angle  should  be  used  for  quantitatively  examin¬ 
ing  inelastic  transition  probabilities:  i t  was  found  that  the  intensities  of  the 
inelastic  peaks  did  not  scale  linearly  with  the  specular  beam  at  two  two  slightly 
different  incident  angles  due  to  the  presence  of  a  selective  adsorption  resonance 
in  the  incident  channel.  Such  resonances  are  common  and  pronounced  features  on 
all  of  the  rare  gas  systems  we  have  examined  [7].  (Our  scattering  geometry  does 
not  eliminate  intensity  variations  that  are  due  to  resonances  in  the  exit  channel.) 

Figure  3  is  a  reduced  zone  plot  of  the  principal  modes  observed  on  the  1,  2, 

3  and  25  layer  layer  Kr  films.  Only  the  energy  loss  and  gain  features  adjacent 
to  the  elastic  TOF  in  Fig.  2  are  plotted.  The  layer-by-layer  evolution  of  the 
surface  dynamics  of  these  structures  is  clearly  revealed,  with  striking  differ¬ 
ences  apparent  between  the  monolayer  and  25  layer  film. 

The  monolayer  SP^mode,  the  only  mode  we  observe,  is  dispersionless  within 
our  experimental  error.  This  indicates  that  the  physi sorbed  Kr  atoms  are  behav¬ 
ing  as  independent  Einstein  oscillators,  vibrating  in  an  uncorrelated  manner  in 
the  well  of  the  holding  potential,  V(z).  We  note  that  this  behavior  has  also 
been  directly  observed  [3,  8,  9],  and  indirectly  inferred  [10],  for  rare  gases 
physisorbed  on  graphite  and  other  metallic  surfaces.  The  monolayer  SPt  frequency 

allowed  us  to  determine  the  second  derivative  of  V(z)  at  the  equilibrium  position, 

.  -2 

and  hence  the  force  constant,  3  171  meV-A  .  Using  this  parameter  in  con¬ 

junction  with  the  previously  determined  well  depth  (D  =  108  mev)  [11],  03  con- 

,3  . 

stant  and  reference  plane  (C^  3  2263  meV-A  ,  zQ  =  0.2521  A)  [12],  we  have  derived 
an  exponential -3  parameterization  of  the  holding  potential  [6],  V(z)  = 

A  exp(-s(z-z  ),-C,/(z-z  )3,  with  A  =  21556  meV  and  g  =  2.407  A'1  (this  potential 
is  not  valid  for  small  z  as  no  saturation  term  was  included  to  account  for  the  sin¬ 
gularity  at  z  ).  This  analysis  demonstrates  that  information  on  phys isorption 


potentials  for  heavy  adsorbates  can  be  obtained  using 
inelastic  He  scattering,  complementing  the  selective  adsorption  techniques 
that  are  applicable  for  light  adparticles  [13].  Fitting  the  bottom  of  the 
well  with  a  Morse  potential  shows  V( z)  to  be  nearly  harmonic  up  to  level 
n  *  5.  Figure  4  is  a  plot  of  all  the  inelastic  peaks  for  the  monolayer.  The 
average  energies  are,  from  bottom  to  top,  2.93,  5.82,  and  8.78  meV— corres¬ 
ponding  to  1,  2,  and  3  times  the  fundamental  within  the  precision  of  our 
measurements .  Since  the  potential  is  harmonic  near  its  minimum,  and  the 
oscillators  are  at  or  near  the  ground  state  for  the  surface  temperature  of 
these  experiments  (25  ±  1  K),  we  can  not  yet  unambiguously  determine  the 
relative  contribution  from  multi  phonon  and  overtone  excitations.  Detailed 
comparisons  between  the  experimental  transition  probabilities  (including 
temperature  dependent  data)  and  those  from  scattering  calculations,  assuming 
either  overtone  or  multiphonon  excitations,  should  help  clarify  this  situation. 

As  further  overlayers  are  added  there  is  increasing  dispersion  (Fig.  3), 
until  at  25  layers  a  well  developed  Rayleigh  mode  is  seen.  The  energy  at  the 
M  point  (short  wavelength  limit)  is  very  nearly  the  bulk  value  after  adsorbing 
only  two  layers.  Not  surprisingly,  the  convergence  to  bulk-like  behavior  is  much 
slower  at  r,  where  displacements  can  extend  well  below  the  surface  [1].  We 
were  unable  to  extend  the  25  layer  measurements  to  Q/Q  <  0.25  because  the 

***  ***1110  X 

weak  inelastic  TOF  peak  merged  with  the  tail  of  the  strong  elastic  signal. 

Also  note  that  not  all  of  the  points  plotted  in  Figs.  3  and  4  are  for  normal 
processes,  G  =  (.00).  Inelastic  transitions  involving  Umklapp  phonons  are  in 
fact  quite  strong  due  to  the  high  corrugation  of  the  He-Kr(lll)  potential. 

All  four  possible  combinations  of  AQ  and  AE  are  observed. 

The  dispersion  curves  were  primarily  mapped  out  by  tracking  the  most 
intense  energy  loss  feature  across  the  first  SBZ  for  Q  <  0  (parallel  momentum 
loss  with  respect  to  the  scattered  He).  The  experimental  intensity  for  this 
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peak  was  ^10  of  the  specular  beam  and  'x-lO"4  of  the  incident  beam,  with  a 
variation  of  no  more  than  a  factor  of  5  across  the  SBZ.  This  small  change 
in  transition  probability  as  a  function  of  Q  differs  from  the  large  expo¬ 
nential  fall -off  found  on  metals  [14],  and  is  partly  due  to  the  larger  "cut¬ 
off  factor",  Qc  [15],  for  this  system  [6,16]. 

Theoretical  dispersion  curves  for  these  overlayers  have  recently  been 
derived  [17,  18,  19].  These  lattice  dynamics  results  fit  the  experimental 
data  quite  well  when  realistic  rare  gas  pair  potentials  and  physisorption 
potentials  are  used  in  the  calculation.  The  agreement  between  theory  and 
experiment  gives  further  evidence  that  the  examined  rare  gas  films  order 
in  a  layer-by-layer  fashion  on  Ag(lll),  i.e.  exhibit  Type  1  growth  [20]. 

We  are  now  examining  deconvoluted  inelastic  transition  probabilities 
and  lineshapes.  In  our  preliminary  analysis  one  interesting  trend  has  already 
been  ohserved  which  is  related  to  phonon  lifetimes:  at  small  Q,  when  the 
adsorbate’s  Einstein  mode  overlaps  the  surface  vibrational  modes  of  the  sub¬ 
strate,  the  linewidths  are  on  the  order  of  1  meV  FWHM;  whereas  at  large  Q, 
where  there  is  no  overlap  and  the  decay  pathway  is  anharmonic,  linewidths 
are  relatively  narrow,  <0.25  meV.  These  observations,  and  the  others  reported 
above,  convincingly  demonstrate  that  inelastic  He  scattering  can  be  used  to 
study  the  surface  dynamics  of  overlayers,  complementing  elastic  diffraction  and 
selective  adsorption  measurements  which  yield  information  on  static  properties 
[7,  13]. 
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Figure  1. 


Figuj'®  2. 


Figure  3. 


Figure  4. 


He  diffraction  spectra  for  Kr  overlayers  physisorbed  on  AgClll), 
taken  along  the  <112>  direction,  with  9.  *  45°,  a  18  meV, 
and  Ts  *  25  K.  The  four  spectra  were  collected  using  the  same 
conditions,  so  the  intensities  are  directly  comparable.  Arrows 
indicate  the  angles  at  which  the  TOF  spectra  of  Fig.  2  were  taken 

Typical  TOF  spectra  for  Kr  overlayers  physisorbed  on  AgClll), 
taken  with  the  same  parameters  as  in  Fig.  1.  Crosses  are  experi¬ 
mental  data,  solid  lines  are  least-squares  fits,  and  dashed  lines 
are  simulations  (.see  text).  Arrows  indicate  the  position  of 
elastic  time-of-f light. 

Phonon  dispersion  curves  of  the  Kr  overlayers  from  r+M.  Open 
symbols  are  for  energy  loss  of  the  beam,  filled-in  symbols  are 
for  energy  gain.  Representative  error  bars  are  shown  for  energy 
and  momentum  loss  features  (with  respect  to  the  incident  beam) 
which  occurred  in  the  first  zone,  i.e.  for  normal  processes. 

Solid  lines  are  polynomial  fits  to  the  data.  The  inset  is  a 
typical  TPD  spectrum  for  a  tri layer.  Such  spectra  were  essential 
as  their  integrated  areas  were  used  to  confirm  the  Kr  coverage  of 
each  film. 

Dispersion  curves  for  monolayer  Kr/Ag(lll),  showing  the  three 
observed  frequencies.  Open  and  filled-in  symbols  have  the  same 
meaning  as  in  Fig.  3. 
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